home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / zlacon.z / zlacon
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. ZZZZLLLLAAAACCCCOOOONNNN((((3333FFFF))))                                                          ZZZZLLLLAAAACCCCOOOONNNN((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      ZLACON - estimate the 1-norm of a square, complex matrix A
  10.  
  11. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  12.      SUBROUTINE ZLACON( N, V, X, EST, KASE )
  13.  
  14.          INTEGER        KASE, N
  15.  
  16.          DOUBLE         PRECISION EST
  17.  
  18.          COMPLEX*16     V( N ), X( N )
  19.  
  20. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  21.      ZLACON estimates the 1-norm of a square, complex matrix A.  Reverse
  22.      communication is used for evaluating matrix-vector products.
  23.  
  24.  
  25. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  26.      N      (input) INTEGER
  27.             The order of the matrix.  N >= 1.
  28.  
  29.      V      (workspace) COMPLEX*16 array, dimension (N)
  30.             On the final return, V = A*W,  where  EST = norm(V)/norm(W) (W is
  31.             not returned).
  32.  
  33.      X      (input/output) COMPLEX*16 array, dimension (N)
  34.             On an intermediate return, X should be overwritten by A * X,   if
  35.             KASE=1, A' * X,  if KASE=2, where A' is the conjugate transpose of
  36.             A, and ZLACON must be re-called with all the other parameters
  37.             unchanged.
  38.  
  39.      EST    (output) DOUBLE PRECISION
  40.             An estimate (a lower bound) for norm(A).
  41.  
  42.      KASE   (input/output) INTEGER
  43.             On the initial call to ZLACON, KASE should be 0.  On an
  44.             intermediate return, KASE will be 1 or 2, indicating whether X
  45.             should be overwritten by A * X  or A' * X.  On the final return
  46.             from ZLACON, KASE will again be 0.
  47.  
  48. FFFFUUUURRRRTTTTHHHHEEEERRRR DDDDEEEETTTTAAAAIIIILLLLSSSS
  49.      Contributed by Nick Higham, University of Manchester.
  50.      Originally named CONEST, dated March 16, 1988.
  51.  
  52.      Reference: N.J. Higham, "FORTRAN codes for estimating the one-norm of a
  53.      real or complex matrix, with applications to condition estimation", ACM
  54.      Trans. Math. Soft., vol. 14, no. 4, pp. 381-396, December 1988.
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.